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calculated. We present the calculation of the anomalous dimensions for all production channels 
of squarks and gluinos and provide numerical predictions for the Tevatron and the LHC. 
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1. Introduction 

Extensions of the Standard Model (SM) involving supersymmetry (SUSY) are one of the best 
motivated models of new physics. The searches for signals of supersymmetry are currently carried 
out by the Tevatron experiments ^ and in the coming years will be undertaken by the experi- 
ments at the Large Hadron Collider (LHC). Within the Minimal Supersymmetric Standard Model 
(MSSM) [|3j], the dominant production processes of supersymmetric particles (sparticles) at hadron 
colliders are those involving pairs of coloured sparticles, i.e. squarks and gluinos, in the final state. 
Depending on the outcome of the experimental searches, predictions for the total rates for these 
production processes are either used to draw the exclusion limits for the mass parameters or will 
help to determine the masses of the sparticles. It is therefore important to know the cross sections 
for hadroproduction of squarks and gluinos with high theoretical accuracy. 

There are four pair-production processes of squarks and gluinos in hadronic collisions: pp — > 
qq ,pp — )• qq ,pp — > gg ,p { p — > qg . The next-to-leading order (NLO) SUSY-QCD corrections [Q, 
H]] to hadroproduction processes are generally positive at the central scale. Depending on the pro- 
cess and the masses of sparticles considered, they can be large. As pointed out in []|], an important 
part of the contributions to the hadronic cross sections comes from the energy region near the par- 
tonic production threshold, reached when the partonic center-of-mass (c.o.m.) energy squared s 
approaches 4m 2 , where m is the average mass of the produced particles. In this region, two types 
of corrections dominate: the Coulomb corrections, due to exchange of gluons between slowly mov- 
ing massive particles, and the soft-gluon corrections, due to emission of low-energetic gluons off 
the coloured initial and final states. The large size of the soft-gluon emission contributions can 
be traced back, for the perturbative n-th order correction, to the logarithmic terms of the form 
a"log /: (j8 2 ) where k = 2n,... ,0 with /3 = \J\ —4m 2 /s. The effects of the soft-gluon emission are 
taken into account to all orders in perturbation theory by performing resummation of the threshold 
logarithms. 

Here we review calculations of the resummed total cross sections for the squark and gluino 
hadroproduction presented in [|| §]. The LO Coulomb corrections to qq and gg production have 
been resummed in [Qj. For the squark-antisquark production process the dominant contribution 
to the next-to-next-to-leading order (NNLO) correction coming from the resummed cross section 
at next-to-next-to-leading-logarithmic (NNLL) level has been studied in Recently, a general 
formalism allowing for simultaneous resummation of threshold and Coulomb corrections has been 
presented and applied to qq production [|T(]] . 

2. NLL threshold resummation 

The inclusive hadroproduction cross section Gh l h 2 ^ki for two massive SUSY particles k and /, 
where k, I can be a squark (q), antisquark (q) or gluino (g), can be written in terms of its partonic 
version Oij^ki as 

Oh^^kiipAtn 2 }) = L / dxidx 2 dp 5 ( p--^— )///ft 1 (^i,M 2 )/;7/ l2 (^2,A t2 )cJ ! 7^/(p,{m 2 },/i 2 ) , 

i j J \ Xl*2 J 

(2.1) 
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where {m 2 } denotes all masses entering the calculations, i,j are the initial parton flavours, f^ hl 
and f jjh 2 the parton distribution functions, and \i is the common factorization and renormalization 
scale. For the hadronic c.o.m. energy y/S, we define the threshold variable p = (m& + m/) 2 / '5 = 
m 2 /S. It measures the distance from threshold for inclusive production of two final-state particles 
with masses m& and m/ in terms of the energy fraction. The partonic equivalent of this threshold 
variable is defined as p = p/{x\X2), where x\^ are the momentum fractions of the partons. The 
resummation of threshold logarithms is performed in the space of Mellin moments N of the cross 



section ( fLl| ), taken with respect to the variable p. Since the soft radiation carries colour charge, 
it can change the colour state of the underlying hard scattering. Consequently, the components 
in the resummed formula describing hard-scattering as well as the components corresponding to 
soft emission carry dependence on the colour configuration of all four particles participating in the 
reaction ij — > kl. This dependence can be vastly simplified if one performs the calculation in an 
orthogonal basis in the colour space for which the soft anomalous dimension matrices, governing 
the soft emission, are diagonal. In this case, up to next-to-leading-logarithmic (NLL) accuracy, the 
resummed partonic cross section reads fllTI |l2| , [l3 , 14 , [TJ|, 16], 



off kl (N,{m 2 },^ 2 ) = [dp p N - 1 a, 7 ^(p,{m 2 },M 2 ) 



:£^4w/(^(M 2 ))A^+l,GV (2-2) 
/ 



The functions A,(iV+ l,Q 2 ,n 2 ) resumming the collinear or collinear and soft radiation from 
the incoming partons are universal and their form can be found e.g. in The new elements which 
have to be calculated in order to obtain the resummed predictions from (Oh are the LO partonic 



cross-sections in Af-space, ofj\ kl j(N, a s (jU 2 )), and the soft function Af"^ kl j(N + 1, Q 2 ,]X 2 ). Both 
types of contributions are dependent on the colour configuration of the underlying hard-scattering 
ij — > kl and have to be calculated separately for each colour exchange channel, indicated by the 
index /. We perform the calculations in the s-channel orthogonal colour basis. 

The soft function, describing large-angle soft gluon emission, at the NLL accuracy is given by 

lnA ij4ki.AN = J Q dz —rr z ^ Dij^ij. 

The coefficients Djj^kij are related to the threshold limit jS — > of the one-loop anomalous di- 
mension matrices Tij^. In this limit, the matrices calculated in the ^-channel orthogonal colour 
basis become diagonal. Each diagonal element Tij^klji contributes to the corresponding Ay^/W,/ 
coefficient. The expressions for m the case kl = qq can be found in [14], kl = gg in [§, 0] 



and kl = qq,qg in []8||. The corresponding expressions for the LO A^-space colour contributions and 
the Djj^kij coefficients have been calculated in [^, [7]] and in [Jsj] for the qq, gg and qq,qg produc- 
tion, respectively. We list the Dij^kij coefficients for all the subprocesses ij — > kl of interest in 
Table 1 . The values of the coefficients equal the negative values of the quadratic Casimir operators 
belonging to the irreducible SU(3) representations for the outgoing state in a given process. This 
confirms the physical picture of the soft gluon radiation at production threshold being governed by 
the total colour charge of the heavy-particle pair in the final state flT0| , |T6| ] . 
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qq^qqj 
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{0,-3} / = {!,»} 

{0,-3} / = {1,8} 
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gg^qqj 



{0,-3,-3} 7 = {1,8 S ,8 A } 

{0,-3,-3,-6,-8} / = {l,8 s ,8 A ,10eT0,27} 

{-4/3,-10/3} / = {3,6} 

{-4/3,-10/3,-16/3} I = {3,6,15} 



D 



D 



gg^ggJ 
,(1) 



D 



qq^qqj 
,(1) 



qg^qgJ 



Table 1: The soft coefficients Dij^ki.i together with the corresponding irreducible SU(3) representations 
spanned by the vectors of the colour basis. 

3. Numerical results 

The resummation-improved cross sections, C7nll+nlo, are obtained through matching the 
NLL-resummed results with the complete NLO cross sections, in the manner described e.g. in 
Ref [J7], |p. The NLO cross sections are evaluated using the PROS PI NO code [fijl], based on 
calculations [|], §]. We assume the left-handed and right-handed squarks in the final state to be 
mass-degenerate and sum over chiralities. No top squarks in the final state are considered. In our 
calculations we use the MSTW 2008 NLO pdfs Jl9| ] and two-loop running strong coupling constant 
a s in the MS scheme with five active flavours. 

We find that the NLL corrections to hadroproduction of squarks and gluinos are positive and 
lead to higher predictions for the total cross sections compared to the NLO results. In general, the 
most significant effects arise for processes characterised by high contributions from gluon initial 
states and the presence of gluinos in the final states. This is to be expected due to the high colour 
charge and, correspondingly, the large values of the Casimir invariants involved. The impact of the 
NLL resummation on the total cross section for the inclusive squark and gluino production, pp — >• 
qtj,qcj,qg,gg + X, is illustrated in Fig. 0. The relative K-factor # NLL - 1 = Cnll+nlo/ Cnlo - 1, 
shown for the Tevatron in Fig. [I]a) and for the LHC in Fig. [j]b), is calculated for various values 
of a fixed ratio r of squark and gluino mass, r = mg/mq. For the average mass m = 600 GeV the 
correction to the inclusive cross section at the Tevatron due to NLL resummation can be as high 
as 18%. The inclusive corrections are smaller at the LHC for sparticle masses below 3 TeV (see 
Fig. [j]b). Given the sparticle mass ranges that we consider, this is consistent with the fact that the 
distance from threshold, i.e. the value of the variable 1 — p = 1 — 4m 2 /S, is on average larger at the 
LHC than at the Tevatron. 

In Figs. |2|a) and^) we show for the Tevatron and LHC, respectively, the resummed NLL+NLO 
total cross section for inclusive squark and gluino production as a function of the average sparticle 
mass m. For illustration we show these results for the choice niq = nig. The error bands indicate the 
theoretical uncertainty of the NLL+NLO total cross section due to the scale variation in the range 
m/2 < \i < 2m. Both at the Tevatron and at the LHC, soft-gluon resummation leads to a reduction 
by around 30% in this part of the theoretical uncertainty. The results presented in Fig. § are the 
most accurate theoretical predictions currently available for the above processes. 
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Figure 1: The relative NLL /^-factor Xnll — 1 = Onll+nlo/onlo ~~ 1 f° r the inclusive squark and gluino 
pair-production cross section, pp/pp — > qq + qq + qg + gg +X, at the Tevatron (a) and the LHC (b) as a 
function of the average sparticle mass m. Shown are results for various mass ratios r = m^jmq. 
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Figure 2: The NLL+NLO cross section for inclusive squark and gluino pair-production, pp/pp qq + 
qq + qg + gg + X, at the Tevatron (a) and the LHC (b) as a function of the average sparticle mass m. Shown 
are results for the mass ratio r = trig /mg = 1 . The error band corresponds to a variation of the common 
renormalization and factorization scale in the range m/2 < fx < 2m. 
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